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D The quark-gluon plasma (QGP) excitations, corresponding to the scalar and pseudoscalar meson 

quantum numbers, for different temperatures are calculated. Analysis is performed in the Hard 

^ Thermal Loop (HTL) Approximation and leads to a better understanding of the excitations of 

QGP in the deconfined phase and is also of relevance for lattice studies. 
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I. INTRODUCTION 



In this paper we consider the scalar and pseudoscalar mesonic correlation function at 
high temperature Quantum Chromodynamics (QCD) in the framework of the Hard Thermal 
Loop (HTL) approach. The evaluation of the mesonic correlator at finite momenta allows, 
by Fourier transform, to get information on its large distance behavior. This, in turn, is 
governed by the mesonic screening mass, a quantity which has been evaluated in lattice 
QCD at large temperature. Hence the present, perturbative results can be compared and 
tested versus the lattice results. 

A similar approach was carried out by the authors a few years ago for the pseudoscalar 
channel only[T]. The results were encouraging but showed some non-negligible discrepancy 
with respect to the lattice data. The present work contains significant improvements in 
the numerical precision of the calculations as well as a new approach to the analytical 
representation of the spatial correlators, thus providing a more accurate extraction of the 
mesonic screening masses in the QGP. Moreover the investigation is extended to the scalar 
channel, for which more recent lattice QCD data are available, hence allowing for a more 
complete discussion of the present result. 

The features of the meson-like excitations inside the hot Quark-Gluon Plasma (QGP) 
provide interesting information on the persistence of interaction effects up to rather large 
temperatures, a characteristic which renders QGP a somewhat peculiar status of matter and 
strengthens the interest on the present experiments which aim to investigate and clarify the 
many issues and questions posed by QCD calculations. 

In spite of the well-grounded reputation of lattice QCD results, an analytical calculation, 
like the one presented in this work, although contained within the limits of the model 
approach, may allow a deeper understanding of the physical behavior of quarks inside the 
plasma and the identification of the relevant degrees of freedom. 

The paper is organized as follows: in Section 2 we shortly recall the details of the cal- 
culation of the mesonic spatial correlation function and spectral function, both within the 
HTL approximation and for the free case. In Section 3 we illustrate the details of the fitting 
procedure adopted in order to derive precise values of the asymptotic mass. The results 
for the scalar and pseudoscalar channels are shown and compared with a few lattice data. 
Comments and conclusions are reported in Section 4. 



2 



II. SCALAR MESON SPATIAL CORRELATION FUNCTION 



The spatial correlation function is conveniently obtained from the finite temperature 
correlator of currents carrying the proper quantum numbers to create/destroy mesons: 

Gm{-it, x) = {Jui-ir, x)jliO, 0)) , (1) 

where r G [0, (3 = 1/T] and Jm denotes the fluctuation of the current operator Jm = q TmQ, 
the vertex Tm selecting the appropriate channel (scalar, pseudoscalar, etc.). The correlator 
in Eq.Q is usually expressed through its Fourier components: 

ImHt,x) = ^ E y J^,e-'^''^e'P■'^GM{^cOn,p) , (2) 



G J 



Un = 2mTT {n = 0, ±1, ±2 . . . ) being the bosonic Matsubara frequencies. 
We are now interested in the 2-axis correlation function G{z), 



Gm{z) = j '^^ J dx±GMi-ir,x±,z) , (3) 



where the integrations select in Fourier space the vanishing frequency and transverse mo- 
mentum components. The asymptotic behavior of static plane-like perturbations with the 
quantum numbers of a (scalar) meson is expected to be exponentially suppressed at large 
values of z (the distance from the plane in the transverse direction) 

g{z) ~ e-"'--^ . (4) 

z—^+oo 

The suppression parameter, which is governing the large distance behavior, is the so-called 
screening mass rriscr- It gives information on the nature of the excitations characterizing 
the QGP phase; the correlator g{z) is frequently studied in lattice simulations, from which 
the screening mass can be extracted as well. 

One can express the spatial correlation function through the finite momentum spectral 
function of the quark- ant iquark excitations in QGP, a{uj,p). Here we shall consider explicit 
formulas for the scalar channel only, though results will be presented both for the scalar 
and pseudoscalar channels. The latter was dealt with in Refs [H |2] and we refer the reader 
to these works; moreover in the massless case chiral symmetry provides similar formulas for 
the scalar and pseudoscalar cases, as it will be pointed out below. 



In the scalar channel the spectral function as{ou,p) is related to Gs{z) by: 



+00 



-00 
+00 



where Pz is the momentum in the z-direction. This is the starting point for investigating 
2;- axis correlations. 

In spectral representation one can relate the spectral function as{u},p) with the meson 
propagator xs{i'^n,p) as follows: 

+00 

Xsii(^n,p) = - [ du^^:^^—^ as{uj,p) = -Imxsiuj + ir],p). (6) 

J lOJn TV 



-00 



The scalar meson 2-point function p], in turn, is 

Xs{ii^hP) = -NfN,- 22 / j^Tr[l S{iuJn,k)l S{iuJn- iuJi,k - p)] , (7) 

n=— 00 ^ ' 

where S{iUn, k) is the quark propagator and NfNc are the numbers of flavors and colors. 

In HTL approximation the quark propagator is dressed by the interaction with the other 
particles of the thermal bath (antiquarks and gluons), carrying typical hard momenta (pro- 
portional to the temperature of the thermal bath); in the spectral representation one can 
write: 

+00 

S {ibjn,p) = - / du— , (8) 

—00 

where the HTL quark spectral function is 

Phtl(w,P) = ^ P+{(^,P) + ^ , (9) 

with 

P±(^, k) = ^^^-^IH^ - + + + 0±(^, k)e(k^ - . (10) 

= -^77^ (11) 



and 



2 [k{-u±k)+ml{±l-^ln'^)Y+[-^ml^] 
Here the thermal gap mass of the quark is niq = gT / a/6. 



Inserting ([s]) into ^ and setting q = k — p one gets: 



hoo +00 



1 



-00 —00 



X Tr[lpHTL(wi,fc)lpHTL(w2,g)] 



(12) 



Then we sum over the Matsubara frequencies in Eq. ( 12 ) with a standard contour integration, 



we perform the usual analytical continuation iuji — )■ w + irj^ (corresponding to retarded 
boundary conditions) and extract the imaginary part of the result, thus obtaining: 



roo +00 



—00 —00 

X S{uj - ui - UJ2) ■ Tr[lpHTL(wi, A;)1phtl(-c^2, q)] ■ (13) 



Now, by inserting Eq. ^ into Eq. (13) and since 

Tr 



Tr 



[l-k.q) 
[1 + k.q) 



(14a) 
(14b) 



one gets 



d^k 
(2^ 



+00 



(e^'^ — 1) / dui I du2n{uji)n{u2)S{uj —cui —UJ2) 



—00 —00 



X |(1 + fc • q)[p+{uJi,k)p+{uJ2,q) + p- (wi, A;)p-_ (^2, g)] 

+ {l-k- q)[p+{ui,k)p-{u2,q)+ p-{u:uk)p+{u2,q)]^ , (15) 

which is exactly equal to the o"ps(w,p) of Ref. [2], since for massless quarks chiral symmetry 
holds (the thermal mass mq acquired in the bath does not affect the symmetries of the 
original Lagrangian). 

In the case of a non-interacting system of quarks with mass m, at finite momentum the 
free quark spectral function reads: 



pf^'^K) = {If+m)^[5{k, - efc) - 5{k, + e,)] 



(16) 
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where = a/Zc^ + and K = (ko,k). Then, by inserting [instead of p^"^^] (16) into 



Eq. (13) one gets the analytic expression 



free i \ 



c^V/, ,2 2 



87r2 



— p — 4m 



X < 6'(w -p - 4m 



4m2 2 , 



^2 _ p2 pP 



+ (^ip'-^\pB\, (17) 



where 



A = log 1 



-loK 1 + e 



= log 1 + e 



-log 1 + e 



For high energy (w — ?■ oo) the spectral function diverges quadratically (cr — )■ w^), hence in 
order to get a finite result for Q{z) Eq. (|5]) one has to regularize integrals. For the non 
interacting QGP the problem was solved for all mesonic channels in Ref. [3] by adopting the 
Pauli-Villars regularization scheme. In the scalar channel one gets: 
NfN.T 



r> f ree 



+00 



(18) 



X 



^ -2z^{2l + + m^ (^2zv/(2/ + + m^ + l) 



l=—oo 



The same procedure cannot be adopted for the interacting case, since the meson spectral 
function is obtained only via numerical evaluation. Nevertheless a "numerical" regularization 
of the integrals can be performed by following the same method proposed in Ref. pLj. We 
notice that the high frequency divergence is of the same order both for the free and for the 
interacting system; hence we first define the difference between the two spectral functions: 



a 



diffi 



The numerical evaluation of this quantity was accurately tested to converge to zero for high 



energies and, for a given momentum p^, it can be numerically integrated as well. In Eq. (19) 



the asymptotic mass m = \/2Tng has been used in the free spectral function (the so-called 
auxiliary spectral function of Ref. [1] ) . 

Here we shall not limit ourselves to consider the scalar channel, for which no result exists 
yet: we also present results for the pseudoscalar channel since we have significantly increased 
the precision of the numerical calculations with respect to previous works [1]. 
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FIG. 1: The meson scalar and pseudoscalar differences of spectral functions, Eq. (19) divided by u, 
for different momenta p = Pz and different temperatures: T/Tc = 1 in the upper panel, T/Tc = 4 
in the lower panel. 
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In Fig. [T]we show the u dependence of Eq. (19) for a few vahies of the momenta pz- The 
numerical calculations were carried out up to a; = in order to check the perfect 

convergence to zero, in the infinite uj limit, of the difference between the full and the free 
spectral functions. 

Then we fit the numerically obtained Qg^^^{pz) = Qs^'^^{Pz) — Qs{Pz) to a sum of Yukawa 
like functions: 

e*"(p.) = i:JT^5- (20) 

i=l '""i ^ ^2 

which can be easily converted by the Inverse Fourier Transform (IFT) to C^"^^^ {z) in coor- 
dinate space 

+ 00 

gd^ff{^^)= I ^ ^tPzZ gd^ff^p^^^J2c,e-"'^^ . (21) 

— oo 

The parameters rrii and q for different temperatures ranging from 1 to 10 Tc are collected 
in Table |T] for the pseudo-scalar channel and in Table |TT] for the scalar one. We remark that 
two terms are sufficient in the pseudoscalar channel, while the scalar one requires up to four 
terms, depending upon the temperature. The number of digits reported in the tables is 
required to ensure an accurate fit of the function. 



T/Tc 


Cl mi 


C2 m2 


1 


-0.7159 6.7027 


-23.5261 148.9711 


2 


-3.5432 13.5179 


-106.3161 285.0631 


4 


-17.7817 24.6371 


-290.4082 251.5351 


10 


-202.8622 60.8533 


-1931.3612 355.4231 



TABLE 1: Parameters of Q'^'^-ff (p^) defined in Eq. (20) for the pseudoscalar channel; the mi are in 



The accuracy of the fit is shown in Fig. [2] for the two channels, at two different tem- 
peratures above Tc. One can notice a similar behavior (apart from sign and size) at the 
higher temperature, while T/Tc = 2 shows that the two channels have the same sign, but 
the minimum in the scalar channel is absent in the pseudoscalar one. 



After Fourier transforming, according to Eq. (21), the final formula for the full spatial 
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T/T, 
1 
2 
4 
10 



Cl 



C2 



m2 



C3 



C4 m4 



-92.31 84.98 -834.71 512.86 -17.92 17.52 -1.92 5.34 
-44.95 15.86 30.08 8.42 -323.78 317.84 
104.97 10.27 1088.71 487.25 
1176.17 231.28 1084.73 30.03 1081.82 30.05 88.52 10.54 



TABLE II: Parameters of G'^^^^iPz) defined in Eq. (20) for the scalar channel; the rrii are in [fm 



correlation function in coordinate space is: 



(22) 



For large z the final results ( 22 ) is once again fitted to an analytic function as follows 



(23) 



i=l 



At variance with the previous approach, where we have assumed that the interaction 
part G'^^^-^Ipz) of the correlation function is dominating and we have extracted the screening 
mass TUscr = "tTT-i dircctly from its asymptotic contribution, namely from the lightest masses 
of Table |l] (only the pseudoscalar channel was considered in Ref. [1]), here we improve the 
fitting procedure by extracting the screening masses directly from the fit of the the total 



spatial correlation function Eq. (23). This makes a small but significant difference (in the 



previously considered channel): indeed now the masses approach the high temperature limit 
from below the free, massless non interacting system limit {2itT). We should also notice 
that, while in the pseudoscalar channel the choice of the lightest mass is rather obvious at 
all considered temperatures, in the scalar channel this procedure would be questionable (see 
Table Inl). 

In Table |V] we summarize the results for the screening masses at all temperatures and 
for all channels considered here. The numbers in the table are calculated assuming Nf = 2 
and Tc = 202 MeV for the transition temperature, according to Ref. [6] (this value is also 
close to the one resulting from the lattice calculations we shall compare later on). The 
ratio Tc/Ajjg = 0.7721 for the Nf = 2 case was taken from Refs. \T-'W\. Throughout the 
calculation we have employed a running gauge coupling as given by the two-loop perturbative 
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FIG. 2: Fit for the meson scalar and pseudoscalar difference, Eq.(20) of the momentum p 
correlation function as a function of w, for two values of T/Tc- 



T/T, 


bi mi 


b2 


m2 


1 


6.843 6.953 


2.381 10~2 


4.875 


2 


26.643 13.505 


8.133 10-3 


10.466 


4 


101.705 26.389 


-5.991 10-1 


23.027 


10 


626.345 65.553 







TABLE III: Parameters of Q^^^\z) defined in Eq. (20) for the pseudoscalar channel; the rrij are 
[fm-i]. 
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FIG. 3: Fit for the meson scalar and pseudoscalar (23) spatial correlation function for different 
T/Tc values. 



beta-function, leading to the expression: 



g''{T) = 2bo log + log {2 log -^|, 
Where bo = d^{n- 2^), 61 = ^{l02 - 38^) 



(24) 



(16,^1 ^ iju^ The choice of the renormalization 
scale // should reflect the typical momentum exchanged by the particles which, in an 
ultra-relativistic plasma, is of order T. Here, for the sake of simplicity, we adopt the choice 
H = 1.1427rT, which was suggested in Ref. [B]. For reference, we also recall that the thermal 
gap mass of the quark is = g{T)T/\/Q. 



Next we compare our results for the screening masses with recent lattice results. In 
particular we refer to the lattice data of Cheng et al. [10], which are the most recent and 
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T/T, 


bi mi 


62 m2 


1 


5.8132 6.9986 


0.8941 4.5283 


2 


-18.2976 9.5012 




4 


-66.1941 21.7911 




10 


-166.7731 75.5921 328.4512 64.2514 



TABLE IV: Parameters of Qf'^^\z) defined in Eq. (20) for the scalar channel; the rrij are in [fm 



T/T, 


m/2nT 




mfy2'KT 


mf jlTxT 


1 


0.229 


1.101 


0.758 


0.704 


2 


0.176 


1.060 


0.814 


0.738 


4 


0.150 


1.044 


0.895 


0.847 


10 


0.129 


1.033 


1.019 


0.998 



TABLE V: The columns display, respectively: m = y^m^, the HTL asymptotic quark mass, 
micr^ = 2-v/7r^T^ + m^, the free screening mass, mp^ and m^^^ , which are the screening masses of 
the interacting QGP, in the pseudo-scalar and scalar channels, respectively 



employ two different lattice sizes, but are obtained in pseudoscalar channel only; in addition 
we show a comparison with the results of Mukherjee [11], who also provides data in the scalar 
channel. Fig. |4] shows the temperature dependence of the screening masses obtained in the 
present work, both in the scalar (black circles) and pseudoscalar (black diamonds) channels, 
divided by the free massless asymptotic limit {2nT). Three values of temperature are shown 
(1, 2 and 4 Tc), the last being the upper limit for which lattice data are also available. We 
report our result at T = Tc only for sake of comparison, but we remind the reader that the 
HTL approach can be considered as a reliable one only at higher temperatures (e.g. T = 2Tc 
or more). 

In the left panel we compare with the results of Cheng et al. [TU] in the pseudoscalar 
channel. The agreement is not very satisfactory and appears to be better for the smaller lat- 
tice (open squares). In the right panel we compare with the results of Mukherjee [H], which 
allow for a full comparison in both channels: with the exception of the lowest temperature 
points, already commented above, the agreement between lattice data and our calculation 
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is rather good. The scalar and pseudoscalar lattice results are closer to each other than in 
our approach, but the trend and the size are quite comparable. 
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FIG. 4: Temperature dependence of the scalar and pseudoscalar screening masses compared with 
the lattice results extracted from Refs [TUl E] 

The same considerations apply to the results illustrated in Fig. [5} where we extend 
the scale of the temperature up to 10 Tc, the highest value considered in this work: here 
we also compare with the screening mass of the "non-interacting quarks" , namely with the 
screening mass value of the auxiliary spectral function employed to evaluate the full HTL 
results: mll^^ = 2v^7r2TM-^T?. The latter is obviously always larger than 2ttT but we can 
assume it as the asymptotic value, which in the figure appears to be almost reached at 
T = 10 Tc. 



III. CONCLUSIONS 



The main goal of this paper was to calculate the asymptotic mass, at different tempera- 
tures, of the scalar and pseudo-scalar mesons in an interacting quark gluon plasma described 
within the HTL approximation. For this purpose we considered the scalar and pseudo-scalar 
mesonic correlation function at high temperature QCD. The evaluation of the mesonic cor- 
relators allowed us to obtain information on its large distance behavior, which, in turn, is 
governed by the mesonic screening mass. 

The main issue of the present approach, which extends to the scalar channel the method- 
ology employed in a previous work pLj, was to improve the precision in the numerical deter- 
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FIG. 5: Temperature dependence of the scalar and pseudoscalar screening masses compared with 
the lattice results extracted from Ref. [llj and with the screening mass in the non-interacting 
system, m{cr'^ = 2\/7r^T^~+~m^. 

mination of the correlators at large distances, hence obtaining a more precise determination 
of the screening mass, the parameter which governs the large distance behavior. We have 
thus modified the procedure for fitting the numerically determined HTL correlator: the full 
(interacting) correlator is determined via the asymptotic (analytic) behavior of the free sys- 
tem, which allows for an appropriate regularization of the otherwise divergent integrations. 
Then a sequence of fitting procedures, first in momentum space and then in coordinate 
space, allowed us to determine the screening mass more accurately than in the past. 

Although the differences do not appear to be dramatic, an important feature of the 
screening masses presently obtained is that they approach, with increasing temperature, the 
non interacting limit from below rather than from above, as in the previous work. The 
comparison with the available lattice data confirms this characteristic. Moreover the values 
for the screening masses we obtained here are in fair agreement with at least one set of 
lattice data [H]. Only close to the critical temperature our approach shows disagreement 
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with respect to lattice results, which does not come as a surprise, since the HTL approach 
is considered to be reliable only well above T^. 

We also notice that other perturbative computations of the static correlation lengths [T^ 
[T3] , led to a small but positive correction to the free screening mass: 



A being an additional correction term, to the same order. The numerical values turn out 
to be above free limit 27rT , whereas all available results from lattice collaborations lie 
below it [HHHl HOj [TTl [T11423] . This fact strengthens the validity of the present approach and 
particularly of the double fitting procedure which allowed us to obtain smaller screening 
masses within, basically, the same HTL approach used in the past. 
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